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Abstract 

We employ an effective gravitational stellar final collapse model 
which contains the relevant physics involved in this complex phenom- 
ena: spherical radical infall in the Schwarzschild metric of the homoge- 
neous core of an advanced star, giant magnetic dipole moment, mag- 
netohydrodynamic material response and realistic equations of state 
(EOS). The electromagnetic pulse is computed both for medium size 
cores undergoing hydrodynamic bounce and large size cores under- 
going black hole formation. We clearly show that there must exist 
two classes of neutron stars, separated by maximum allowable masses: 
those that collapsed as solitary stars (dynamical mass limit) and those 
that collapsed in binary systems allowing mass accretion (static neu- 
tron star mass). Our results show that the electromagnetic pulse spec- 
trum associated with black hole formation is a universal signature, 
independent of the nuclear EOS. Our results also predict that there 
must exist black holes whose masses are less than the static neutron 
star stability limit. 
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1 Introduction 



What is the measurable signal for gravitational collapse? It certainly pro- 
duces gravitational waves, but unfortunately these have shown to be unde- 
tectable so far. However, there is another signal. Since the star has a sizable 
magnetic field, and since we expect this field to get quenched very rapidly 
during gravitational collapse, a massive electromagnetic pulse (MEMP) is 
created. We further show that the power spectrum is unlike any other nat- 
urally occurring phenomena: a wave packet whose wavetrain lasts only 2 
milliseconds and whose characteristic power spectrum is a square block. 

We present a simulation, within an effective model, of the iron-core col- 
lapse of massive stars which produce type II supernova. Depending on the 
core's mass, either a hydrodynamic bounce occurs above nuclear densities 
or the collapse races unhindered into the formulation of a stellar size black 
hole. Either case produces a MEMP. Our main contribution in this paper 
is the calculation of the energy spectrum of the MEMP. In cases where the 
collapsing core has mass greater than the dynamical neutron star mass limit 
(discussed herein), a hydrodynamic bounce cannot occur and no visible su- 
pernova results. In such cases, the MEMP is the only electromagnetic signal 
associated with the gravitational collapse. 

One of the unexpected secondary results is the finding that, in general, 
there exist two classes of neutron stars: the first class is made up of solitary 
neutron stars which do not have the opportunity to accrete matter from a 
binary companion, and therefore have original masses coming from the col- 
lapse process. This mass is bounded by the ability of the nuclear equation 
of state to produce a hydrodynamic bounce. An unexpected secondary find- 
ing is that this imposes an heretofore unappreciated constraint on candidate 
nuclear equations of state. The maximum mass of the neutron star in this 
class is the dynamical mass, being the largest core mass that can undergo 
a hydrodynamic bounce. We will see that this dynamical mass limit is con- 
siderable less than the static mass limit, so the bounce does indeed save the 
core from further collapse due to its own self-gravity. The second class of 
neutron stars are those that can accrete matter from a binary companion 
and so their maximum mass is constrained by the requirement of static sta- 
bility: the maximum static neutron star mass. In general, we find that the 
maximum static mass can exceed the dynamic mass by as large as 100%. 

A current review of the status of gravitational collapse is given in Joshi 
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(2000), who discusses spherically symmetric collapse. In the present paper, 
we consider spherical radial infall. As shown in the work of Bocquet et 
al. (1995), only gigantic magnetic fields B > 10 10 T = 10 GT (testla = 
10 4 gauss, GT = 10 9 T) cause stellar deformation. Since the magnetic fields 
considered in this paper fall far below this large critical field, the star is 
not deformed by its magnetic field. An interesting followup to the research 
presented here, is to apply the effective gravitational collapse model to stars 
having non-negligible rotation. 

In general, there are three sources for electromagnetic and gravitational 
radiation associated with stellar collapse: the direct radiation phase emit- 
ted by the stellar object before the formation of a black hole, the so-called 
damped oscillations known as quasi-normal ringing (Iyer 1987) that are the 
vibrations of a black hole shaking off its non-zero multipole charge moments, 
and the late-time power-law tail (Price 1972; Cunningham 1978; Leaver 1986; 
Hod 1999, 2000) underneath the damped oscillations. Because of the gravi- 
tational redshift, essentially the only energy that survives to infinity comes 
from the direct radiation phase, before the black hole forms. In this paper, 
we calculate the direct phase electromagnetic energy radiated by stellar ob- 
jects that bounce and become stable neutron stars, and stellar objects so 
massive that they become black holes. 

In order to calculate the electromagnetic radiation produced by the grav- 
itational collapse, our strategy is the following. During the collapse the infall 
kinetic energy decreases while the system does mechanical work against the 
internal star pressure. This pressure has, in general, two components: mate- 
rial (nuclear) pressure and electromagnetic pressure. However, all published 
nuclear equations of state have the property that the nuclear pressure is or- 
ders of magnitude larger than the electromagnetic pressure; indeed, only for 
magnetic fields greater than 10 18 Gauss does the electromagnetic pressure 
approach parity with the nuclear pressure. Using conservation of energy, 
the continuity equation and the material equation of state, we can calculate 
how the star radius changes with time. If we further assume that the stellar 
magnetic field can be approximated by a dipole field, and that the lines of 
magnetic field are frozen in the material and carried along with it, we can 
relate the change in the star's radius with time to the change in its effective 
magnetic moment. The usual dipole radiation formula then gives the energy 
spectrum of the radiation. 
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Figure 1: Difference between proper and Schwarzschild times. 



We simulate the collapse to the hydrodynamic bounce, or to twice the 
gravitational radius, a g = 2GM/c 2 , if no bounce occurs. In the former case, 
the bounce cuts off the MEMP, while in the latter case, final black hole forma- 
tion at radius = a g occurs almost instantaneously later after the simulation 
ends. Because of the gravitational redshift, the final stellar trajectory from 
2a g to a g adds little to the observable energy at infinity. Thus the MEMP 
calculated is expected to be an accurate profile. The reason why we stop at 
2a g is because, while the Schwarzschild metric is the exterior metric during 
gravitational collapse and as such, is the clock time in the effective collapse 
model, the radiation process must be calculated in proper time. In Fig. 1, 
we plot the difference between proper time and Schwarzschild time. Only 
when the core radius becomes less than 3a g is there a difference in the two 
clocks. By stopping the simulation at 2a g , we can still reliably calculate the 
core surface as a function of proper time, by using the Taylor expansion em- 
ploying the Schwarzschild time solution of the collapse model. Indeed, one 
can calculate from this figure that the time error is 8.364 microseconds x X, 
where Ai/Ai Q = X, with M. being the core mass. At the same time, we fold 
in the gravitational and Doppler redshift, since the difference between the 
Schwarzschild time and proper time leads to a gravitational redshift of the 
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produced radiation. 

The collapse process from the initial starting configuration of the iron-core 
mass density (p) of p = 1 x 10 12 gm/cm 3 and inward speed v/c ~ 2/3 x 1CT 2 
to either hydrodynamic bounce or twice a g takes less than 3 milliseconds. In 
contrast, the characteristic time scale for neutrino diffusion out of the core is 
~ 2 — > 10 seconds. This 10 3 times larger neutrino diffusion scale means that 
these particles cannot influence the bounce or the MEMP. The role of the 
neutrinos is not to carry away the binding energy of the collapsing matter, 
as it sometimes is erroneously stated. Rather, the collapsing matter must do 
work against the equation of state and this energy is stored as internal en- 
ergy, entirely analogous to a weight compressing a spring. This stored energy 
is not lost, but remains primarily as elevated Fermi levels of electrons and 
nucleons. Thus the energy the neutrinos actually carry away is the energy 
from nuclear beta decay, relaxed elastic energy and thermal heat energy. The 
difference in the time scales for neutrino diffusion and hydrodynamic infall 
means that the iron-core collapse is adiabatic, allowing a concise formulation 
for the conservation of energy. Previous work, Morley (1999), showed that 
neutrino transport in the iron-core was due to elastic scattering and typically 
~ 1000 collisions occurred before escape. Our collapse code in the present 
paper shows beyond a doubt, that if no hydrodynamic bounce occurs, these 
diffusing neutrinos are trapped in the infalling nuclear material and are car- 
ried into the resulting black hole. Thus observation of neutrinos associated 
with particular supernova proves that the initial infall did not form a black 
hole. 

Another physical insight into the problem is the recognition that the 
iron-cores have flat density profiles. In Fig. 2, we show the neutron star den- 
sity profile using the Friedman-Pandharipande equation of state, Straumann 
(1992). Thus the approximation that the star collapses homogeneously is 
excellent (t = Schwarzschild time): p(r,t) = p(t). In this paper, we treat 
spherical radial infall. As mentioned earlier, a very interesting follow-on 
to this research is to remove this constraint and allow core rotation in the 
effective collapse model. 

Further simplifications can be achieved by using the fact that the electrical 
conductivity of the iron-core is very high, nearly infinite. This allows the 
magnetohydrodynamic approximation: the lines of magnetic force are frozen 
into the material and are carried along with it. In this paper, we approximate 
the iron-core stellar magnetic field by a magnetic dipole of moment rh. We 
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Figure 2: Neutron star density profile using the Friedman-Pandharipande 
equation of state, Straumann (1992). 

use the model of a uniformly magnetized sphere of magnetization per unit 
volume M. The relation of m to M is 

m = ^-a 3 M (1) 

— * — * 

where a is the stellar radius. B out = H outside the sphere and is a magnetic 
dipole field: 

^ 47T ,, „-*., r 2 cos 6* „ sinfl-., , , 

B out = -jAM\{— —r + —9} (2) 

while inside the core 

B m = yM . (3) 

This is shown in Fig. 3. In general, a = a(t) and M = M(t). The 
magnetohydrodynamic material response is the constancy of the magnetic 
flux through the polar region, Fig. 4 

7ra 2 (^) sin 2 0B in (ti) = 7ra 2 (t 2 ) sin 2 6B m (t 2 ), \B in \ = B in (4) 

giving 

a 2 (ti)|M(ti)| = a 2 {t 2 )\M{t 2 )\ = constant = C m . (5) 
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Figure 3: Stellar magnetic dipole field. 



Therefore 




(6) 



and 



rh(t) = — C m d(t), m = \rh 



(7) 



C m is fixed by astrophysical data. Young pulsars (e.g. PSR B0154+61) have 
the polar field strength B = 2GT. Since this is over a flux area of 10 x 10 
km 2 , the constant is C m ~ 8 x 10 25 cm 2 -Gauss. 

Thus we have obtained a relation between the second order time derivative 
of the stellar effective magnetic moment and the star's surface acceleration, 
which will allow us to obtain the power radiated once we know the accelera- 
tion. This will be obtained by a numerical simulation, whose ingredients we 
explain in the next sections. 

2.1 Conservation of Matter 

The conservation of baryon number is given by the equation 



dr 



= —n y -u 



(8) 
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Figure 4: Constancy of the magnetic flux through the polar region. 

where ribar is the baryon density, r is proper time, y is the 4-derivative and 
u is the 4-velocity. As long as r ~ t and extreme relativistic speeds are 
not encountered, then the conservation of baryon number reduces to the the 
conservation of matter: 



The only region in the collapse phenomena where the two conservation laws 
begin to differ is near twice a g . Though we treat p as uniform throughout 
the spherical volume, it is changing in time t so % = & 7^ 0. Thus v cannot 



V-W + -f = 0. 



(9) 



be a constant. With V • (pv) = pV • v = 
conservation of matter equation becomes 



(r 2 v r ) where v = v r r, the 




(10) 



This has the physical solution 



r dp 
3p dt 



(11) 
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2.2 Conservation of Energy 



The energy available to do work if the matter collapses from mass density 
Pi to p f is 

Eif = PEi + KEi — PEf (12) 

where KE iy KEf is the kinetic energy at the initial, final density respectively, 
and PEi, PEf is the gravitational self-energy of the core with initial, final 
mass density respectively. The work (W) > / done against the equation of 
state in compressing the matter is 

r v f 

W = - PdV, P = pressure. (13) 

J Vi 

Here V i: V) are the initial and final core volumes respectively. The core will 
bounce when or if it reaches a density such that 

E if = W. (14) 

For realistic equations of state, we expect a bounce for core masses of the 
order of a solar mass. 



2.3 Gravitational Self-Energy 

A core of radius a in the Schwarzschild metric has gravitational potential 
energy, Weinberg(1972) 

PE= [ a 4Trr 2 {l , 1 }p e dr (15) 

JO 1 A _ 2GM(r) ^ e ^ > 



where c is speed of light, p e is the energy density of matter and Ai(r) is the 
amount of mass within radius r. Since M.{r) = ^M. and p e = 3 ^ c 3 , this 
reduces to 

PE = p e [ a 4nr 2 {l }dr (16) 

which can be integrated to give 

PE = Mc 2 + ^^Wl -Ka 2 -^-= arcsin (17) 

K CL\/K 

where k = ^§£r, G = gravitational constant. 
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2.4 Kinetic Energy of Infall 

The equation of state (EOS) of matter (to be discussed in more detail later) 
can describe "hot" or "cold" nuclear matter. In general, the thermal energies 
that are reached in the collapse phenomena are a small fraction of the grav- 
itational binding energy experienced by the particles. If the pressure P in 
Eq. (13) describes a "hot" EOS, then some of the kinetic energy + potential 
energy of infall is turned into heat. The numerical procedure (to be discussed 
in detail in a later section) takes small time slices (At) and computes the 
new kinetic energy after work against the EOS is determined. The equation 
is 



f V(t+At) 

KE(t + At) = KE(t) + PE(t) - PE{t + At) + ^ PdV . (18) 



Thus the effective gravitational model takes into account any desired 
heat generation through the inputed EOS. We need the instantaneous kinetic 
energy of a collapsing iron-core in the Schwarzschild metric. We want the 
expression to be accurate to 0(v 4 /c 4 ). 

To order 0(v 2 /c 2 ), the instantaneous kinetic energy is 

r 2/Pe{r)y r {r) 1 

I\ Hi = 7 <±7( T 

Jo 

where 



KE= / ATtr\'^^) r \> , =dr (19) 

V c 2 r 



B = O 2 ■ (20) 



This becomes 

B f a r 4 dr 



r< ra r dr 

= 7T^' (21) 

Computing this integral: 

Q fl 3 3 3 

KE = 27rp e — { y/l — kcl 2 -aVl — K,a 2 H ^-j- arcsin a^/H} . (22) 



c 



4k 8k 2 8k 5 / 2 



The kinetic energy to 0(-u 4 /c 4 ) is 



KE 



= r^^m - , 1 dr . (23) 

Jo K 2 J[ C 2 8 C 4 J / 1 2G^(r) v ; 
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This extra piece can be integrated in the same manner as the first piece, 
and using the parameters a (speed of surface) and A = na 2 we find that the 
kinetic energy to 0(t> 4 /c 4 ) is 

KE = ^{^(^^[-^A^v^-^v^A + ^arcsinv/A] 

+ Y ( ^ [ ^ ( M + ^ + T20 ) - 64 arCSmv/X] (24) 
2.5 Nuclear Equation of State (EOS) 

The hydrodynamic bounce is independent of the equation of state below nu- 
clear density (nuclear density ~2x 10 14 gm/cm 3 ). We know this because if 
we put in zero pressure in Eq. (13) below nuclear density, the bounce prop- 
erties hardly change at all. The hot thermal gas kinetic EOS below nuclear 
density is not responsible for the hydrodynamic bounce, and this can be 
shown by calculating the work W. In this paper, we use two completely dif- 
ferent nuclear EOS, and while the contribution of the pressure below nuclear 
densities is negligible, we still hook them to EOS of degenerate electronic 
matter used frequently in collapse calculations, Baym (1971). 

The different regimes for the nuclear EOS are listed in Table 1. 



Class 


Property 


Applicability 


normal matter 
iron-core 
neutronic 


N ~ Z 
N ~ 2Z 
N » Z 


atomic nuclei 
proto-neutron star 
cold neutron star 



Table 1: The three regimes for the nuclear EOS. The first and third are in 
beta decay equilibrium, but the middle is out of equilibrium. Here N = the 
number of neutrons and Z = the number of protons. 

In this paper, two different nuclear EOS were used in the collapse calcula- 
tions: the "Argonne AV14 + UVII", Wiringa (1995) and the QCD, Kislinger 
(1978). 

Perturbation theory in quantum chromodynamics (QCD) breaks down at 
~ twice nuclear density, where it must be joined to phenomenological EOS. 
Both EOS used here are joined to the Baym EOS, as mentioned earlier. At 
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gure 5: Code logic. 



4 a 
V 



dp* 
~dt. 



a d 2 p 



(25) 
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3 Radiation 



The MEMP calculated here is magnetic dipole radiation; it is not the emission 
of a black-body at some temperature and radius. The MEMP, while it exists, 
completely swamps thermal emission. Our goal is to calculate the MEMP 
energy spectrum. The energy radiated in all directions per unit time, in 
Gaussian units, is, Landau(1975): 

f = 3>MI 2 • < 26 > 

We introduce the Fourier transforms m(uj) 

m(r) = -L= / m(u)e iuJT duj (27) 

m{u) = f m(r)e- iuJT dr . (28) 

Then 

m(r) = ]= f u 2 m{u)e iWT du . (29) 

Taking the integration over time from — oo to +oo gives the total energy 
radiated 

E = ^J u 4 \m{uj)\ 2 duj (30) 

and the energy spectrum 

d l = ^ W «)l' ■ (3D 

3.1 Roll-on-Roll-off (RORO) Function 

The core's dipole moment is 

m(r) = yC m a(r) (32) 

where a(r) is the core's radius in km at proper time r. As discussed in the 
later sections, we fit a(r) to a polynomial curve and compute m(r). Then we 
integrate analytically to obtain the total energy release. To know the energy 
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spectrum, though, and to include the Doppler and gravitational redshifts, 
we must compute the Fourier transform. If we put in Heaviside functions 
in time for the beginning and ending of the infall, the Fourier transform 
will not give rise to a finite energy spectrum because the derivative of a 
step function is a delta function. Physically we must impose the following 
boundary conditions: from time — oo we had the initial radius and at time 
+00 we have the final radius, and smooth transition between them in the 
infall period. The analogous problem occurs in particle scattering. Thus we 
need a Roll-on-Roll-off (RORO) function to accomplish this mathematically. 
Any suitable RORO function will do and gives the same answer. Thus the 
time dependent m(r) becomes (with Co = -f-Cm) 

1 r 
m(r) =>• C -a max (1 - tanh(— )) 

T T - T 

+.25(1 + tanh(— )) x (1 - tanh( ))m(r) 

/IS fJS 

1 r-T 

+C -a min (l + tanh( )) . (33) 

In the above equation, m(r) on the right-hand-side is the fitted poly- 
nomial function in proper time r (fifth order polynomial - fitted from the 
collapse code output of the acceleration), r and T = collapse time, are in 
microseconds, and a max and a min are respectively the starting and ending 
core radius in km. Any derivation of the gravitational or electromagnetic 
power spectrum must use a RORO function. As we will see in the following 
sections, the RORO function is also critical to deriving the redshifts. 

3.2 Total Energy 

We fit the acceleration collapse solution to a polynomial function 

5 

a(T)=km^ft(»)W ( 34 ) 

i=0 

where r is in microseconds. Then 

5 

m( T ) = C km J2 b(i) (t//is)' (35) 

i=0 
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Figure 6: Integration contour, 
must compute m{uj) where m{uj) is given by Eq. (28). Letting 
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uj • /is, x = t//is and Q = T/ fis, we have 



m(uj) = us km C ^= / dx e~ lsx {. ^-IUffLn _ tanhx) 
V27T J-oo km 

5 

+.25(1 + tanhx)(l - tanh(a; - Q)) £ 6(i)x* 



+ 5 ^ (l+tanh(a ,_g ))} (3g) 

This integral can be done using contour integration in the complex energy 
plane. Letting x = Z r + iZi (complex), we have damping for Zi < 0. Poles 
in the lower half of the complex plane occur at Z r = 0, Zi = — 2A ^ +1 ir and at 
Z r = Q, Zi = — 2jV 2 +1 7r. We integrate clockwise along the contour of Fig. 6. 
to obtain the answer: 

/iskmCo , U A _ S7r 2iv+i 
mM = - — 7^—(-2m){}^ e S7T * [-.5(a max /km) 



2vr N=0 

1 + tanhg^ . , 2iV+l 

+ ' 25 tanhQ gfeO)(-^(^-)) J ] 

+ g e - fa Q- W r(^ )[ _ 25 l+tanhQ * 

(Q - irf(^^)y + .5(a mm /km)]} (39) 

The presence of the exponentials signals that the integration over u in Eq. 
(30) is indeed finite. 



3.4 Redshift 

The difference between the proper time and the Schwarzschild time only 
occurs when the radius is below three times a g . In order to compute the 
gravitational and Doppler redshifts, we need to separate out the contribution 
of each accelerated surface to the radiation spectrum. This is done by the 
use of the RORO function. Consider the first RORO function from r = 
to t = Ti, and the second RORO function from r = T\ to r = r 2 . One 
can show that the complex T\ poles of these two RORO functions mutually 
cancel in their sum, so the sum of two RORO reduces to one RORO from 
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Figure 7: Radii time trajectory, for core masses of 1.5 and 5 solar masses, 
for the initial conditions stated in the text. Also, an acceleration profile of a 
core bounce in the Argonne EOS, for a 0.75 solar mass core. 

In the collapse studies, whenever a star goes below three times a g , we 
solve for the radius as a function of proper time by means of the Taylor 
expansion 

a(T) = a(t) + f t (T-t) + .5^(T-tf (40) 
So a(r) is reconstructed from three times a g to either the bounce radius or 
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Figure 8: Spectrum of the radiated energy, for the Argonne (for a core of 
1.0825 solar masses) and QCD (for a core of 1.3806 solar masses) equations of 
state, and for a core mass of 1.5 and 5 solar masses, for the initial conditions 
mentioned in the text. 

It is also of interest to see the kinetics of black hole formation. Unex- 
pectedly, we find that once the maximum core mass allowing hydrodynamic 
bounce is exceeded, the black hole is formed independently of the nuclear 
EOS. Thus, for example, the radius function a(r), infall speed v r (a), and ac- 
celeration function acc(r) are all identical for a 1.5 solar mass core collapse, 
whatever nuclear EOS is used. Thus the MEMP of black hole formation is 
a universal signature. Though this has been derived in an effective gravita- 
tional collapse model, we believe this is a true theorem in physics. It is of 
interest then, whether the introduction of stellar rotation also reintroduces 
nuclear EOS dependence for black hole MEMP. It is probably true that once 
the core mass exceeds the neutron star dynamical mass, whether it is rotating 
or not, the black hole MEMP is a universal function. In Fig. 7, we give the 
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kinetics of a 1.5 solar mass core collapse and a 5.0 solar mass core collapse, 
along with an example of a core bounce. 



4 Results 



Core 


bounce ? 


Initial radius 


Final density 


Final radius 


Time 


0.75 


yes 


70.88 


6.40(14) 


8.21 


2175 


1.00 


yes 


78.01 


7.80(14) 


8.48 


2145 


1.38 


barely 


86.87 


7.64(14) 


9.50 


2110 



Table 2: Results for the QCD equation of state. Core is in M , Radius in 
km, time in /xs, density in gm/cm 3 . 6.40(14) means 6.40 xlO 14 . The static 
neutron star mass for QCD is 2.34 solar masses. 



Core 


Emitted MEMP Energy 


Received MEMP Energy 


QCD 0.75 


1.07(43) 


1.07(43) 


QCD 1.00 


1.47(43) 


1.47(43) 


QCD 1.3806 


2.05(43) 


1.98(43) 


ARG 0.75 


1.27(43) 


1.27(43) 


ARG 1.0825 


1.82(43) 


1.76(43) 



Table 3: The emitted and received (difference due to redshifts) directly ra- 
diated MEMP energy of the QCD and Argonne (ARG) EOS. Units are in 
ergs. 

The code logic is given in Fig. [5|. Following Arnett (1977), the iron core 
begins collapse with a fairly uniform density of 10 12 gm/cm 3 and initial infall 
speed v/c ~ | x 10~ 2 . After choosing the mass of the core, the initial radius 
is determined. 

In Fig. 7, we present the results we obtained for the stellar surface for core 
masses of 1.5 and 5.0 solar masses, and in Fig. ^ the results for the energy 
spectrum of the emitted electromagnetic radiation, for the Argonne EOS (for 
a core of 1.0825 solar masses) and QCD (for a core of 1.3806 solar masses) 
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Core 


bounce ? 


Initial rad. 


Final dens. 


Final rad. 


Time 


0.75 


yes 


70.88 


1.60(15) 


6.05 


2195 


1.0825 


barely 


80.10 


1.29(15) 


7.36 


2144 



Table 4: Results for the Argonne EOS. The static neutron star mass for 
Argonne is 2.10 solar masses. 



Core 


Initial rad. 


Time 


Rad energy 


1.5 

5 


89.3 
133.40 


2142 
2001 


2.66 (43) 
7.60(43) 



Table 5: 1.5 and 5 core mass collapse 

EOS cases, and for core masses of 1.5 and 5.0 solar masses (independent of 
EOS). 

5 Conclusion 

We have presented an effective gravitational collapse model that is thought 
to include the essential physics: spherical radial infall in the Schwarzschild 
metric of a homogeneous core of an advanced star possessing a giant mag- 
netic dipole moment. The electrodynamic equations are represented by the 
approximation of magnetohydrodynamic material response, which will be an 
excellent approximation due to the high electrical conductivity. Any nuclear 
EOS can be used, and we chose two which have high interest: the "Ar- 
gonne" AV14+UVII which represents the typical phenomenological nuclear 
EOS and the QCD EOS. Surprisingly, the maximum core mass that bounced 
(the dynamic neutron star limit) is rather small using the phenomenological 
EOS. 

The main result of the paper is the calculation of the MEMP from the 
collapse. Both the total energy and its spectral characteristics are derived. 
In cases where the core collapse is greater than the dynamical neutron star 
mass, the MEMP is the only electromagnetic experimental signature. It 
is hoped that this will spur activity to develop a suitable experimental re- 
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ceiverfj. To answer the question as to whether the surrounding material from 
the collapsing object will quench or change an observational signal is clearly 
complicated, and it requires a detailed calculation in transport theory us- 
ing information on the supernova models that vary tremendously on their 
precursor environment, which is beyond the scope of the present paper. 

The existence of two classes of neutron star masses leads to the situation 
where a stellar core greater than the dynamical mass, but less than the static 
mass, collapses to form a black hole, thus producing black holes less massive 
than the mass limit of neutron star stability For example, a black hole of 
1.6 solar masses can exist. This possibility sheds light on one of the most 
perplexing problems of astrophysics: How do black holes of stellar masses 
form? As envisioned here, a core more massive than the dynamical mass 
collapses and forms a black hole. It may now accrete the remaining outlying 
surface matter to produce a whole continuum of stellar mass size black holes. 

After this paper was submitted for publication, we were informed of the 
work of Hanami ||1997| . This author attempted to explain gamma-ray bursts 



using the change in the magnetic field of a collapsing stellar object. Unfor- 
tunately, there are several problems associated with this work. Hanami did 
not consider the fact that the star radiates continuously throughout the col- 
lapse trajectory and his numerical solution violates Maxwell's equations. He 
used a zero pressure EOS, which will not give the correct surface acceleration 
and core remnant. In order to compute the total energy, the gravitational 
redshift must be properly included. The calculation of the gravitational red- 
shift is a very difficult problem as explained in the present paper, because, 
as the star collapses and radiates, each instantaneous stellar configuration 
has a different redshift. In order to do this, one needs to employ multiple 
Roll-on- Roll-off mathematical functions in the complex energy plane, correct 
for the gravitational redshift and add the spectrums incoherently. Hanami 
did not consider this and simply calculated the redshifted power at the ini- 
tial collapse configuration, to obtain his result. At the bottom of page 688 
and top of page 689, Hanami [|1997|| alludes to the fact that his solution is 
the quasi-normal ringing (QNR), but his eigenfrequency is incorrect (QNR 
have precise discrete eigenfrequencies) . Also his time dependence does not 

1 Due to the peak of the spectral curve around the wavelength ~ 2km, the receiver will 
have to be a satellite able to detect a broadband spectrum within the electromagnetic 
pulse. The additional benefit is that terrestrial signals in this frequency regime cannot 
penetrate earth's ionosphere, thereby reducing spurious man-made noise. 
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correspond to any known quasi-normal mode damping. 
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